Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1199 1592 (N36W28 L=293)  1398    358        1 (1037/1537-1340/1622) -12.3 W15-W38 N32-N39
      2 1524 1601 (N27W53 L=317)   854    124        1 (1454/1578-1590/1621) -12.1 W46-W60 N25-N29
      3 1646 1529 (N19W59 L=323)   542     71        1 (1612/1496-1681/1561) -78.7 W55-W61 N16-N21
      4  207 1307 (N39E66 L=198)   497     38        0 ( 185/1289- 235/1324)   7.1 E69-E61 N38-N40
      5 1451 1274 (N11W32 L=296)   269    110        0 (1409/1236-1496/1303) -51.3 W29-W33 N07-N13
      6 1055 1209 (N17W06 L=270)   431    152        0 (1017/1136-1079/1257)  82.9 W04-W07 N12-N20
      7  527 1071 (N20E30 L=234)   327     83        1 ( 499/1050- 568/1103) -38.7 E32-E27 N17-N21
      8  521  921 (N11E33 L=231)   415    162        0 ( 486/ 844- 562/ 991) -73.3 E34-E31 N05-N15
      9 1878  819 (S28W61 L=325)  2906    212        6 (1816/ 746-1933/ 910)  18.4 W52-W73 S32-S24
     10  641  753 (S01E29 L=235)   920    309        0 ( 576/ 638- 746/ 820)  21.8 E36-E21 S06-N02
     11 1578  795 (S20W29 L=294)   234     54        0 (1554/ 771-1604/ 820)  26.6 W27-W31 S20-S18
     12  482  638 (S05E43 L=221)   870    217        2 ( 408/ 593- 570/ 684) -35.0 E48-E37 S09-S01
     13  462  552 (S10E48 L=216)   690     79        3 ( 421/ 520- 499/ 586)   0.0 E52-E45 S12-S08
     14  306  518 (S11E67 L=197)  2016    129        5 ( 256/ 469- 355/ 567) -50.2 E74-E60 S14-S07
     15  589  520 (S14E40 L=224)   365     81        0 ( 566/ 482- 611/ 564)  80.5 E41-E38 S17-S11
     16  519  469 (S17E48 L=216)  1006    216        4 ( 476/ 399- 564/ 546)  36.9 E55-E41 S20-S11
     17 1447  424 (S42W15 L=279)  1160    210        3 (1364/ 358-1532/ 505)  17.4 W07-W23 S45-S36
     18  597  406 (S22E45 L=219)   694     90        1 ( 555/ 362- 622/ 448) -80.5 E47-E41 S25-S19
  

    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