Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1214 1837 (N53E21 L=292)   519     60        0 (1182/1816-1244/1857)  21.8 E24-E17 N50-N55
      2  646 1801 (N30E58 L=255)  1020    132        1 ( 605/1741- 695/1854) -57.5 E62-E54 N24-N35
      3  866 1763 (N34E37 L=276)  1165    213        2 ( 840/1658- 898/1844) -33.7 E46-E28 N26-N39
      4  879 1440 (N13E20 L=293)  1141    302        4 ( 838/1302- 910/1585) -64.8 E26-E16 N04-N21
      5  397 1393 (S01E50 L=263)   307     50        0 ( 373/1370- 421/1413)  33.7 E51-E48 S03-N00
      6 1426  612 (S16W37 L=350)   529     97        1 (1384/ 589-1479/ 646)  38.7 W34-W40 S18-S14
      7 1607  537 (S14W54 L=  7)   388     33        0 (1585/ 512-1628/ 558)  90.0 W52-W55 S15-S12
      8 1240  514 (S28W30 L=343)   920    158        3 (1155/ 496-1327/ 538)  37.9 W25-W35 S30-S23
      9 1497  452 (S22W51 L=  4)   999    130        3 (1428/ 417-1553/ 499)  49.4 W47-W54 S25-S18
     10  516  371 (S65E23 L=290)   394      7        0 ( 496/ 355- 537/ 386)  26.6 E27-E19 S67-S62
     11  727  307 (S61W08 L=321)  2427    320        4 ( 599/ 270- 849/ 358)  28.4 E09-W18 S68-S52
  

    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