Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1234 1669 (N40E08 L= 96)  2437    475        9 (1066/1532-1346/1780) -19.1 E24-W02 N31-N43
      2  638 1606 (N18E42 L= 62)   650    115        3 ( 598/1547- 673/1660) -82.9 E45-E40 N13-N21
      3  429 1473 (N05E51 L= 53)   335     89        1 ( 409/1438- 466/1514)  80.5 E53-E49 N02-N08
      4  871 1340 (N08E17 L= 86)   667    218        1 ( 824/1241- 910/1451) -45.0 E22-E12 N02-N13
      5  950 1404 (N14E15 L= 89)   741    229        0 ( 893/1339-1033/1451)  40.6 E17-E09 N08-N16
      6  222 1238 (S13E61 L= 43)   947    146        0 ( 194/1159- 250/1302)  83.7 E63-E59 S17-S08
      7  167 1264 (S12E69 L= 34)  1281    157        4 ( 130/1204- 196/1301)  71.6 E73-E65 S15-S09
      8 1685 1168 (N22W37 L=140)  1148    158        3 (1643/1117-1722/1236) -40.6 W32-W40 N17-N25
      9  391 1119 (S16E42 L= 62)   457    103        1 ( 354/1078- 425/1164)  66.8 E43-E39 S19-S12
  

    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