Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  389 1385 (N30E48 L=190)  2168    491        2 ( 230/1307- 556/1448)   4.9 E67-E32 N24-N34
      2 1002 1395 (N28W01 L=239)   366    182        0 ( 914/1362-1066/1417) -24.0 E05-W04 N25-N30
      3 1419 1369 (N23W30 L=268)   253     94        0 (1384/1339-1471/1395) -33.7 W27-W33 N21-N25
      4  480 1321 (N25E38 L=201)   627    116        1 ( 418/1300- 545/1338)   0.0 E43-E32 N24-N26
      5  866 1227 (N18E09 L=229)   196    103        1 ( 828/1199- 899/1269)  45.0 E11-E06 N16-N20
      6  205 1174 (N16E64 L=175)  1729    218        3 ( 175/1078- 233/1267) -65.6 E70-E59 N10-N21
      7  337 1172 (N16E49 L=190)   296     78        1 ( 295/1155- 374/1187)   8.1 E52-E45 N15-N17
      8 1327  940 (S02W18 L=257)   695    171        1 (1256/ 881-1367/1003)  40.6 W13-W21 S05-N01
      9 1805  916 (S08W56 L=295)   785    103        2 (1781/ 868-1829/ 967)  63.4 W53-W59 S11-S05
     10 1745  859 (S11W50 L=289)   502     60        1 (1723/ 826-1766/ 891)  76.0 W48-W52 S13-S09
     11  781  559 (S24E20 L=219)   797    179        1 ( 760/ 479- 810/ 635)  65.6 E21-E16 S29-S18
     12  545  550 (S24E38 L=201)   345    107        0 ( 497/ 523- 589/ 582)  20.6 E41-E33 S25-S21
  

    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