Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1323 1833 (N54W36 L=284)  1983    140        2 (1267/1794-1397/1884) -56.3 W29-W43 N49-N60
      2 1518 1573 (N31W39 L=287)  2819    442        4 (1317/1501-1714/1627) -10.6 W22-W56 N26-N35
      3 1677 1367 (N17W47 L=295)   723    159        2 (1648/1298-1712/1442) -83.7 W45-W48 N13-N22
      4 1009 1371 (N18E00 L=248)   631    229        2 ( 921/1305-1107/1433) -35.0 E05-W05 N14-N22
      5 1476 1315 (N14W30 L=278)   414    126        0 (1437/1270-1523/1371) -50.2 W27-W32 N11-N17
      6 1786 1237 (N09W54 L=303)   425     74        4 (1761/1208-1817/1260) -11.3 W51-W58 N07-N10
      7  203  815 (S12E62 L=186)   945    166        3 ( 178/ 737- 234/ 903)  78.7 E63-E60 S16-S06
      8 1720  760 (S20W49 L=298)   417    152        1 (1696/ 687-1756/ 823)  71.6 W47-W51 S24-S15
      9 1779  777 (S18W55 L=303)   432     67        1 (1753/ 746-1807/ 808)  45.0 W52-W57 S20-S16
     10 1314  463 (S39W21 L=269)  1082    228        0 (1220/ 422-1412/ 524)  19.7 W13-W27 S42-S34
     11 1750  480 (S38W68 L=317)   601     68        1 (1728/ 461-1783/ 509)  -7.6 W62-W77 S39-S35
     12 1586  464 (S39W47 L=295)   528     74        1 (1563/ 426-1615/ 508) -40.6 W42-W50 S42-S36
  

    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