Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  589 1540 (N34E37 L=331)   396     95        0 ( 572/1482- 605/1580) -54.5 E39-E33 N29-N36
      2 1584 1480 (N33W43 L= 50)   763    146        2 (1511/1438-1634/1518) -29.1 W36-W46 N30-N36
      3  866 1478 (N31E13 L=354)   621    244        0 ( 779/1404- 945/1514) -21.0 E19-E06 N26-N33
      4 1348 1478 (N33W22 L= 29)   238     85        0 (1313/1446-1379/1506) -38.7 W19-W24 N30-N34
      5 1309 1435 (N30W19 L= 26)   319    206        0 (1217/1389-1386/1464) -18.4 W11-W24 N26-N31
      6 1103 1199 (N14W04 L= 12)   335     79        1 (1073/1177-1146/1225) -11.3 W02-W07 N12-N15
      7  707  967 (S02E20 L=348)   281    131        1 ( 659/ 935- 752/ 987) -18.4 E22-E16 S03-N00
      8  468  934 (S05E36 L=331)   488     70        1 ( 444/ 891- 493/ 971) -78.7 E37-E34 S08-S03
      9 1681  905 (S03W45 L= 52)   803    120        3 (1631/ 867-1735/ 947) -16.7 W40-W50 S05-N00
     10 1048  800 (S11W02 L= 10)   275     94        1 (1009/ 767-1087/ 829) -36.9 W00-W04 S12-S09
     11  435  769 (S16E40 L=328)   289     98        1 ( 418/ 733- 455/ 812)  68.2 E40-E37 S18-S13
     12 1808  769 (S12W61 L= 68)   481     51        2 (1790/ 742-1832/ 798) -53.1 W58-W64 S13-S09
     13 1402  769 (S12W25 L= 33)   251     71        1 (1367/ 747-1439/ 782)  11.3 W22-W27 S13-S10
     14  857  677 (S20E10 L=358)  1015    215        7 ( 798/ 597- 904/ 769) -70.0 E13-E06 S24-S13
     15  965  694 (S18E03 L=  5)   620    110        4 ( 922/ 642-1009/ 739) -50.2 E05-W00 S21-S15
     16  410  568 (S30E46 L=321)  1512    248        3 ( 358/ 508- 508/ 655) -32.5 E49-E37 S34-S24
     17 1138  488 (S32W10 L= 18)  4101    618       13 ( 926/ 389-1374/ 648)  28.9 E04-W25 S40-S19
     18  747  582 (S27E17 L=350)   821    185        3 ( 683/ 533- 823/ 629)  32.5 E22-E11 S30-S23
     19 1497  496 (S30W38 L= 46)  1540    309        4 (1370/ 429-1627/ 553)  25.3 W29-W48 S35-S25
  

    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