Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1225 1468 (N32W11 L=208)  1574    457        2 (1097/1329-1393/1565) -31.8 W01-W23 N23-N38
      2  522 1381 (N21E38 L=159)   569    148        1 ( 498/1324- 569/1446) -45.0 E41-E33 N17-N25
      3 1706 1315 (N24W50 L=247)   461     74        2 (1689/1273-1725/1351) -90.0 W48-W52 N21-N26
      4 1700 1155 (N14W47 L=244)   399     73        1 (1671/1126-1730/1184) -31.0 W44-W49 N12-N15
      5 1867 1122 (N12W66 L=263)   985    114        3 (1842/1063-1890/1179)  54.5 W62-W70 N08-N15
      6  865 1096 (N06E10 L=186)   300    109        0 ( 844/1057- 902/1147) -51.3 E11-E07 N03-N08
      7 1750  857 (S04W53 L=250)   515      8        0 (1723/ 834-1773/ 877)  14.0 W50-W55 S05-S03
      8 1105  699 (S18W07 L=204)   695    244        2 (1019/ 609-1157/ 779)  57.5 W02-W10 S24-S12
      9 1318  643 (S20W22 L=219)   270     84        0 (1287/ 610-1352/ 683)  59.0 W20-W23 S22-S17
     10 1687  607 (S21W54 L=250)   374     34        0 (1670/ 579-1705/ 632) -56.3 W51-W55 S22-S19
     11  367  567 (S33E51 L=146)   687     94        2 ( 320/ 551- 414/ 588)   9.5 E57-E45 S33-S30
     12 1242  415 (S37W22 L=218)   383     71        0 (1207/ 398-1288/ 433)  18.4 W19-W25 S38-S35
  

    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