Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1308 1730 (N54W19 L= 83)   334     22        0 (1284/1717-1332/1743) -11.3 W16-W21 N52-N55
      2  597 1472 (N25E35 L= 29)  1592    357        4 ( 509/1392- 658/1606) -86.2 E41-E30 N19-N34
      3 1501 1402 (N30W31 L= 96)  2467    545        3 (1284/1232-1628/1528) -31.6 W14-W40 N19-N37
      4 1657 1265 (N22W43 L=108)   316     48        0 (1636/1236-1679/1291)  36.9 W40-W45 N20-N23
      5 1365 1191 (N15W20 L= 85)   384    150        1 (1332/1141-1407/1260) -54.5 W17-W23 N11-N18
      6 1544 1196 (N17W33 L= 98)   344    106        0 (1523/1139-1557/1245) -71.6 W31-W34 N13-N19
      7  445  782 (S19E38 L= 27)   693    164        0 ( 422/ 710- 481/ 843)  71.6 E39-E34 S23-S14
      8 1458  650 (S18W33 L= 98)   849    198        2 (1379/ 583-1530/ 708)  48.4 W28-W37 S22-S13
      9 1718  650 (S16W56 L=120)   441     47        0 (1701/ 621-1738/ 675)  71.6 W53-W57 S17-S14
     10 1724  566 (S21W61 L=126)  1744    173        4 (1662/ 514-1799/ 612)  19.7 W53-W72 S25-S17
  

    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