Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  825 1872 (N56E51 L=203)   666     66        0 ( 788/1857- 860/1883)   0.0 E57-E45 N54-N57
      2 1157 1734 (N50E05 L=249)   378     57        0 (1137/1703-1186/1755)  68.2 E06-E02 N47-N52
      3  472 1679 (N32E62 L=192)   960     97        2 ( 442/1646- 505/1722) -11.9 E74-E55 N29-N34
      4 1308 1635 (N45W11 L=265)   426     62        0 (1278/1614-1339/1657) -21.8 W07-W13 N43-N46
      5  769 1563 (N29E28 L=226)   378    124        0 ( 714/1544- 829/1580)  29.7 E31-E23 N27-N31
      6 1794 1355 (N33W58 L=312)   517     45        1 (1761/1338-1818/1378)  12.5 W53-W62 N31-N35
      7 1585 1251 (N23W34 L=288)   675    213        1 (1505/1187-1649/1318) -26.6 W27-W39 N19-N25
      8  292  935 (S17E50 L=204)   350     69        1 ( 267/ 916- 324/ 954)  21.8 E52-E47 S18-S15
      9  275  873 (S21E52 L=202)   507     61        2 ( 253/ 845- 292/ 909) -63.4 E54-E50 S22-S18
     10  715  474 (S41E13 L=241)   256     63        1 ( 684/ 458- 749/ 495)   0.0 E15-E10 S41-S39
     11 1465  439 (S29W45 L=299)  1110    174        1 (1380/ 417-1547/ 471)  17.1 W38-W51 S32-S26
     12 1050  195 (S59W29 L=283)   693     74        1 (1028/ 162-1071/ 236) -33.7 W23-W35 S62-S54
  

    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