Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  856 1860 (N41E48 L=346)  1061     92        3 ( 816/1822- 893/1890)  -8.7 E55-E41 N38-N43
      2  445 1627 (N14E59 L=335)  1987    259        2 ( 393/1516- 495/1751) -50.5 E69-E54 N06-N23
      3 1232 1599 (N34E06 L= 29)   758    134        2 (1206/1537-1258/1667) -52.1 E09-E02 N29-N38
      4  479 1518 (N08E49 L=345)  2120    312        3 ( 368/1426- 600/1613)  56.3 E56-E44 N00-N16
      5  266 1153 (S18E55 L=340)   802    133        0 ( 207/1118- 313/1194)  31.0 E61-E50 S20-S14
      6  122 1013 (S28E74 L=320)   770     58        3 ( 100/ 988- 152/1034)   9.5 E80-E67 S29-S27
      7  587  842 (S28E22 L= 13)   824    201        4 ( 507/ 797- 689/ 901)  42.3 E27-E15 S32-S22
      8  987  632 (S29W10 L= 44)   764    393        1 ( 815/ 602-1168/ 657)  34.4 W00-W19 S35-S22
      9 1398  601 (S18W36 L= 70)   942    197        2 (1316/ 555-1481/ 644)  59.0 W31-W39 S22-S12
     10  997  489 (S37W17 L= 51)  1116    283        1 ( 884/ 408-1102/ 554)  66.4 W12-W20 S46-S30
  

    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