Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1189 1256 (N16W05 L=284)  3876    962        7 ( 886/ 879-1460/1588) -33.1 E20-W29 S01-N31
      2 1729 1561 (N50W52 L=331)   650    153        0 (1693/1549-1761/1577)  16.7 W46-W56 N48-N52
      3 1484 1545 (N42W23 L=302)   434    105        1 (1453/1520-1534/1561)  -8.1 W20-W27 N39-N43
      4  292 1341 (N04E58 L=221)   375     79        3 ( 266/1316- 315/1360)  -9.5 E61-E55 N03-N05
      5 1855 1132 (N22W58 L=337)  1077    245        2 (1797/1083-1912/1165)  -4.1 W51-W66 N19-N23
      6  177 1042 (S15E63 L=216)  1354    287        4 ( 140/ 935- 231/1129)  69.0 E66-E57 S22-S09
      7 1822  900 (N07W59 L=338)   458    162        1 (1775/ 869-1849/ 924)  -5.7 W53-W63 N05-N08
      8  446  815 (S25E35 L=244)   749    184        1 ( 359/ 784- 516/ 846)   8.7 E41-E28 S26-S22
      9  979  778 (S17W02 L=281)  1834    508        3 ( 737/ 728-1232/ 844)  22.5 E13-W16 S21-S08
     10 1424  355 (S34W47 L=326)  5615    845       11 (1135/ 310-1645/ 434)  28.5 W26-W62 S43-S24
     11 1061  356 (S43W20 L=299)   380     28        0 (1032/ 339-1084/ 376)  25.5 W18-W22 S44-S41
  

    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