Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1837 1298 (N16W66 L=121)   738    133        0 (1803/1239-1856/1348) -90.0 W62-W68 N12-N19
      2  508 1276 (N20E35 L= 20)   245     27        0 ( 493/1254- 527/1297) -56.3 E36-E33 N18-N21
      3  318  976 (N02E49 L=  6)   441    119        1 ( 264/ 943- 375/1007) -16.7 E53-E43 N00-N03
      4 1249  914 (S03W13 L= 68)   386    141        0 (1227/ 850-1276/ 986) -76.0 W12-W15 S07-N01
      5 1738  932 (S05W49 L=104)   512     66        1 (1699/ 908-1770/ 964)  15.9 W45-W52 S06-S03
      6  795  916 (S02E14 L= 41)   309     80        1 ( 761/ 878- 829/ 952)  53.1 E16-E12 S04-N00
      7  492  863 (S05E35 L= 20)   300    107        0 ( 464/ 833- 510/ 909)  76.0 E37-E33 S06-S02
      8  258  779 (S11E57 L=358)   883    140        3 ( 229/ 719- 285/ 845)  58.0 E60-E52 S14-S06
      9 1192  799 (S10W10 L= 65)   436    142        0 (1123/ 768-1253/ 830)  20.6 W05-W13 S11-S08
     10  209  718 (S15E66 L=349)  2066    217        5 ( 164/ 629- 258/ 812) -74.7 E71-E60 S20-S09
     11  295  674 (S18E56 L=359)   390     72        0 ( 272/ 643- 318/ 707)  33.7 E59-E52 S19-S15
  

    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