Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1071 1868 (N54E31 L= 27)   515     61        2 (1047/1848-1095/1890)  -9.5 E34-E27 N51-N56
      2  411 1654 (N19E64 L=355)   844    170        1 ( 368/1614- 443/1673)  15.9 E68-E60 N16-N21
      3  942 1575 (N27E21 L= 38)   796    259        3 ( 854/1492-1053/1655) -12.5 E29-E11 N23-N30
      4 1342 1557 (N37W06 L= 65)   438     79        2 (1295/1526-1377/1586)   0.0 W02-W09 N35-N38
      5  327 1495 (N08E62 L=357)  1186    212        4 ( 251/1438- 383/1549)  38.7 E69-E58 N03-N11
      6  485 1296 (S01E40 L= 19)   567    241        1 ( 411/1276- 601/1324)  17.1 E45-E32 S03-N03
      7 1766  774 (N02W56 L=115)   655    115        0 (1739/ 735-1798/ 809) -18.4 W52-W61 N00-N03
  

    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