Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  870 1782 (N60W19 L=297)   312     31        0 ( 851/1759- 887/1800)   9.5 W16-W22 N57-N61
      2 1277 1742 (N41W46 L=324)   457     51        0 (1250/1722-1305/1766) -71.6 W43-W48 N39-N43
      3  681 1692 (N59E06 L=272)   750    143        1 ( 608/1658- 741/1719)   0.0 E15-W01 N57-N60
      4 1127 1588 (N37W24 L=303)  1002    242        3 (1026/1546-1237/1622) -14.0 W17-W33 N33-N39
      5  908 1555 (N41W08 L=287)   306     52        1 ( 873/1534- 941/1571) -21.8 W05-W10 N39-N42
      6  496 1094 (N23E32 L=247)   530    128        1 ( 453/1055- 548/1150) -49.4 E34-E28 N19-N26
      7 1506 1030 (S05W28 L=306)   835    238        0 (1482/ 920-1544/1134)  48.0 W23-W33 S10-N00
      8  482 1000 (N18E34 L=244)   719    249        2 ( 391/ 951- 571/1029) -17.1 E41-E28 N13-N19
      9  460  890 (N12E38 L=241)   380     93        1 ( 440/ 846- 482/ 927)  45.0 E39-E35 N09-N13
     10 1729  702 (S32W40 L=319)  2661    370        4 (1609/ 576-1802/ 852)  32.3 W28-W48 S36-S23
     11  587  551 (S11E39 L=240)   361     68        0 ( 562/ 513- 614/ 586) -80.5 E40-E37 S13-S07
     12  441  482 (S12E54 L=225)   844    169        0 ( 396/ 423- 499/ 531) -71.6 E56-E51 S17-S08
     13  412  398 (S18E63 L=215)  1309    153        7 ( 366/ 341- 460/ 443) -66.8 E67-E59 S22-S13
     14  715  330 (S28E41 L=237)  1416    207        3 ( 646/ 294- 814/ 373) -13.0 E46-E32 S30-S23
     15 1002  243 (S42E28 L=251)  1751    282        3 ( 900/ 169-1122/ 324)  12.1 E42-E14 S46-S38
     16  864  154 (S46E49 L=230)   597     60        1 ( 837/ 134- 889/ 185)  45.0 E54-E43 S49-S43
  

    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