Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  925 1919 (N69W16 L=331)  1312    127        3 ( 877/1870- 965/1956)  13.0 W07-W33 N62-N75
      2 1344 1434 (N20W26 L=342)   618    253        0 (1240/1348-1418/1501) -54.0 W21-W29 N13-N25
      3  843 1480 (N30E06 L=310)   405    111        0 ( 788/1463- 898/1499) -14.0 E09-E01 N28-N31
      4  457 1274 (N22E34 L=281)   407     68        1 ( 438/1238- 480/1318)  63.4 E35-E32 N20-N25
      5 1360 1116 (N01W21 L=336)   269     84        0 (1324/1079-1391/1145)  21.8 W18-W23 N00-N02
      6  586  579 (S21E36 L=279)  1842    356        4 ( 539/ 443- 666/ 753)  76.0 E38-E31 S30-S10
      7 1627  558 (S38W40 L=355)  2368    362        3 (1480/ 470-1768/ 657)  17.1 W27-W53 S42-S32
      8 1042  479 (S34E07 L=308)   284     91        1 ( 993/ 467-1090/ 495)  -8.1 E10-E03 S35-S32
      9  710  393 (S35E34 L=281)   290     83        0 ( 697/ 355- 735/ 422)  33.7 E37-E31 S37-S32
     10  926  281 (S47E23 L=292)   307     79        0 ( 889/ 267- 972/ 295)  -7.1 E26-E18 S48-S45
  

    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