Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  971 1559 (N25E19 L= 69)   401    176        0 ( 907/1522-1011/1622)  69.4 E22-E17 N21-N29
      2 1597 1434 (N37W27 L=115)   388     92        0 (1564/1392-1626/1470) -38.7 W23-W29 N34-N38
      3 1423 1347 (N25W15 L=103)   292     92        0 (1396/1306-1450/1394) -26.6 W11-W17 N23-N27
      4 1687 1304 (N31W36 L=124)   284     71        1 (1667/1280-1706/1338) -33.7 W34-W38 N30-N33
      5 1377 1245 (N17W14 L=103)   659    209        2 (1319/1190-1476/1304)  -4.8 W09-W21 N15-N19
      6 1330 1220 (N14W12 L=100)   215     33        0 (1311/1197-1347/1240) -33.7 W10-W13 N13-N15
      7  145 1102 (S21E71 L= 18)   826    101        1 ( 129/1049- 170/1139)  63.4 E73-E66 S24-S18
      8 1536  737 (S06W39 L=127)   307     64        0 (1497/ 722-1574/ 752)  31.0 W36-W41 S07-S04
      9 1514  622 (S13W42 L=131)   621    146        1 (1456/ 562-1563/ 671)  77.5 W40-W44 S17-S08
     10 1414  424 (S27W47 L=135)   516    167        0 (1332/ 394-1485/ 446)  45.0 W42-W50 S30-S23
     11  899  310 (S52W20 L=109)   291     47        0 ( 885/ 284- 913/ 338) -38.7 W17-W23 S54-S50
     12 1039  255 (S51W36 L=124)   513     82        0 (1008/ 222-1074/ 287) -10.3 W31-W42 S52-S49
  

    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