Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  960 1636 (N35E20 L=179)  2959    640        7 ( 710/1516-1195/1776)  -6.5 E45-W00 N30-N40
      2  386 1594 (N20E61 L=137)  1138    309        1 ( 365/1508- 406/1680) -47.7 E68-E56 N14-N25
      3 1269 1502 (N33W06 L=204)   277     57        1 (1243/1477-1303/1518)  53.1 W03-W07 N30-N34
      4  598 1454 (N15E37 L=162)   380     77        0 ( 575/1422- 632/1499)  80.5 E37-E35 N12-N18
      5  362 1279 (N00E49 L=149)   491     87        1 ( 328/1254- 395/1313)  36.9 E51-E46 S01-N02
      6 1864  961 (N12W62 L=260)  1561    255        3 (1807/ 857-1900/1037) -36.9 W53-W67 N06-N16
      7  740  968 (S11E16 L=182)   232     53        0 ( 713/ 940- 768/ 998) -26.6 E17-E13 S11-S09
      8  901  916 (S11E05 L=193)   404     59        1 ( 860/ 888- 940/ 939)  36.9 E07-E03 S12-S08
      9 1657  686 (S08W50 L=248)   701    170        1 (1640/ 623-1683/ 752) -45.0 W46-W54 S11-S04
     10 1313  633 (S19W27 L=225)   242     37        1 (1286/ 616-1338/ 647)   0.0 W24-W28 S19-S17
     11  987  581 (S29W08 L=206)   409     86        1 ( 941/ 558-1030/ 607)   9.5 W04-W10 S30-S27
  

    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