Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1279 1449 (N30W20 L=299)   561    123        1 (1240/1393-1307/1505)  81.9 W17-W22 N26-N33
      2  460 1380 (N28E41 L=238)   324     90        1 ( 416/1358- 506/1402) -18.4 E45-E36 N26-N29
      3 1367 1355 (N23W25 L=304)   438     66        0 (1330/1335-1410/1378) -26.6 W22-W28 N21-N24
      4  439 1210 (N17E40 L=239)   306    103        1 ( 418/1155- 463/1255)  81.9 E41-E38 N13-N20
      5 1271  830 (S09W15 L=294)   477    106        1 (1243/ 793-1298/ 884)  59.0 W12-W16 S11-S05
      6 1703  804 (S14W47 L=326)  1113    180        2 (1638/ 716-1772/ 883)  39.3 W41-W53 S18-S08
      7 1607  802 (S13W39 L=318)   345    162        0 (1578/ 731-1622/ 867)  77.5 W36-W39 S17-S08
      8 1770  741 (S18W55 L=334)   462     49        1 (1744/ 716-1788/ 762)  14.0 W52-W57 S20-S17
      9 1809  681 (S23W63 L=342)  2195    213        2 (1766/ 587-1854/ 759)  84.8 W59-W68 S29-S18
     10 1731  633 (S25W54 L=333)   735     86        1 (1713/ 587-1754/ 672)  80.5 W50-W56 S28-S22
     11  288  619 (S20E61 L=218)   728     86        1 ( 274/ 574- 304/ 670)  45.0 E64-E56 S23-S16
  

    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