Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1634 1545 (N27W59 L= 46)   454     80        0 (1595/1515-1669/1563) -45.0 W55-W61 N24-N28
      2 1777 1495 (N20W73 L= 61)   844     89        0 (1750/1459-1803/1537) -80.5 W70-W76 N17-N23
      3  719 1395 (N33E16 L=332)   934    224        2 ( 673/1292- 770/1479) -65.2 E19-E12 N26-N39
      4 1094 1381 (N27W10 L=358)   310    119        1 (1056/1332-1139/1426) -60.3 W08-W12 N23-N30
      5 1698 1160 (N03W47 L= 35)   710     93        3 (1660/1120-1726/1201) -78.7 W44-W50 N00-N05
      6 1826 1081 (S04W59 L= 47)   456     85        1 (1811/1035-1837/1119)  68.2 W57-W60 S07-S02
      7  384  697 (S07E49 L=299)   901    151        2 ( 352/ 622- 418/ 761)  41.6 E54-E44 S11-S03
      8 1130  660 (S18W01 L=349)   258    139        0 (1068/ 635-1182/ 694)  15.9 E02-W05 S18-S16
      9 1363  579 (S27W16 L=  4)   235     87        0 (1326/ 562-1407/ 600)   0.0 W13-W19 S28-S26
     10  911  430 (S31E18 L=330)   275    102        1 ( 899/ 378- 924/ 482)  69.4 E18-E15 S34-S26
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y