Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  479 1703 (N69E42 L=112)   803    104        2 ( 428/1675- 519/1730)   0.0 E56-E32 N65-N72
      2 1472 1588 (N26W49 L=203)  2311    383        3 (1341/1557-1646/1621) -25.6 W38-W63 N18-N30
      3 1592 1481 (N16W51 L=205)   514    130        0 (1545/1456-1638/1502) -33.7 W47-W54 N13-N18
      4  486 1267 (N34E30 L=124)   388    107        1 ( 431/1243- 535/1290) -29.1 E35-E26 N31-N36
      5 1119 1260 (N18W12 L=166)   335    176        1 (1046/1223-1199/1281) -41.2 W07-W15 N13-N20
      6  438 1140 (N27E36 L=118)  1116    295        1 ( 399/1048- 526/1274) -69.4 E38-E29 N20-N36
      7  616 1125 (N22E22 L=132)   274     68        1 ( 588/1102- 650/1152) -53.1 E24-E20 N20-N24
      8 1782 1067 (S12W51 L=205)   418     69        1 (1762/1030-1798/1103) -78.7 W49-W52 S14-S09
      9  281  843 (N12E55 L= 99)  1231    283        1 ( 191/ 812- 372/ 913)   2.7 E66-E45 N10-N15
     10  395  769 (N06E47 L=107)   608    133        1 ( 364/ 700- 421/ 821)  81.9 E48-E44 N02-N09
     11  658  706 (S03E30 L=124)   399     61        2 ( 630/ 670- 687/ 740) -63.4 E31-E27 S04-N00
     12  527  538 (S10E45 L=109)   393     59        2 ( 499/ 515- 562/ 562)   0.0 E47-E41 S11-S08
     13  388  462 (S13E62 L= 92)   716     77        0 ( 369/ 416- 405/ 498)  35.5 E65-E57 S15-S10
     14 1155  412 (S34E10 L=145)   355     49        0 (1122/ 384-1176/ 431) -63.4 E11-E07 S35-S31
  

    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