Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  972 1676 (N49W16 L=349)   410     90        0 ( 947/1626- 992/1710)  82.9 W13-W17 N44-N52
      2 1032 1650 (N45W19 L=352)   294     43        0 (1006/1627-1058/1674) -76.0 W17-W20 N43-N47
      3 1195 1610 (N38W28 L=  2)  1204    153        2 (1126/1575-1271/1659) -21.8 W23-W33 N34-N42
      4 1501 1483 (N22W44 L= 17)   649    102        2 (1474/1430-1524/1531)  63.4 W41-W47 N19-N25
      5 1671 1376 (N12W52 L= 25)   556    116        0 (1625/1332-1708/1419)  10.3 W46-W57 N10-N13
      6  462 1325 (N36E34 L=299)  1179    244        2 ( 335/1284- 553/1357) -10.3 E48-E26 N32-N38
      7  730 1241 (N25E14 L=319)   897    238        2 ( 628/1165- 810/1285) -45.0 E20-E09 N19-N29
      8  676 1074 (N16E20 L=313)   370    165        0 ( 637/1030- 730/1102) -31.0 E23-E17 N14-N18
      9  146  942 (N15E71 L=262)  2606    159        4 ( 115/ 866- 198/1020) -63.4 E77-E62 N10-N20
     10  455  935 (N12E38 L=295)   706    186        3 ( 388/ 881- 518/1000) -48.4 E42-E34 N07-N16
     11 1308  967 (S02W15 L=349)   193     46        1 (1293/ 943-1327/ 990)  45.0 W14-W16 S03-N00
     12 1717  895 (S17W43 L= 16)  1046    210        4 (1667/ 839-1789/ 943)  20.0 W38-W50 S20-S12
     13 1239  822 (S09W08 L=342)   351    150        0 (1176/ 783-1288/ 876)  14.0 W04-W12 S10-S06
     14 1125  797 (S08W01 L=335)   190     39        1 (1098/ 776-1149/ 815) -45.0 W00-W02 S09-S06
     15  855  742 (S06E16 L=318)   407    173        1 ( 801/ 709- 911/ 773)  -8.1 E19-E12 S08-S04
     16  884  629 (S14E17 L=317)   212     66        1 ( 843/ 609- 909/ 646)   0.0 E19-E15 S14-S12
     17  587  496 (S17E41 L=293)   317     71        0 ( 554/ 479- 622/ 516) -36.9 E42-E38 S18-S15
     18  634  339 (S29E46 L=287)   388     43        1 ( 604/ 327- 668/ 353) -21.8 E49-E43 S29-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