Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1531 1445 (N14W43 L=143)  1363    473        2 (1335/1371-1730/1523) -42.3 W34-W56 N03-N23
      2 1095 1456 (N27W17 L=117)   380    121        1 (1036/1432-1164/1479)  -6.3 W12-W21 N25-N28
      3 1844 1380 (N00W70 L=170)   882    252        1 (1804/1354-1887/1397)  -6.7 W63-W80 S02-N01
      4 1347 1355 (N14W28 L=128)   340     98        0 (1307/1334-1395/1384) -45.0 W25-W30 N12-N17
      5 1843 1120 (S13W57 L=157)  2317    414        4 (1802/ 944-1895/1269)  72.5 W53-W63 S24-S04
      6  314  958 (N18E47 L= 53)   766    171        4 ( 270/ 906- 369/1003) -49.4 E50-E43 N14-N21
      7 1649  918 (S18W35 L=135)   723    256        3 (1565/ 837-1705/ 986) -65.2 W30-W38 S24-S11
      8 1161  660 (S19E02 L= 98)   388    142        2 (1117/ 620-1211/ 690)   9.5 E05-W01 S20-S16
      9 1294  492 (S34W01 L=101)   623    120        3 (1273/ 441-1318/ 543)  66.8 E01-W03 S37-S29
     10  754  390 (S25E36 L= 64)   923    251        1 ( 717/ 310- 814/ 477)  67.4 E39-E30 S31-S19
     11 1479  389 (S49W13 L=113)  1818    342        7 (1396/ 340-1546/ 454)  -4.1 W05-W20 S54-S43
     12  501  317 (S24E60 L= 40)  1697    314        7 ( 472/ 263- 537/ 378)  38.7 E66-E54 S27-S19
  

    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