Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  827 1737 (N49E08 L= 13)   339     61        2 ( 804/1701- 849/1768) -71.6 E10-E05 N45-N51
      2 1192 1709 (N41W22 L= 43)   355     74        0 (1158/1690-1235/1723)   8.1 W18-W25 N40-N42
      3 1367 1525 (N26W29 L= 49)  1075    229        2 (1249/1473-1449/1551) -26.6 W20-W34 N22-N28
      4 1705 1405 (N16W52 L= 73)  1169    235        4 (1651/1359-1753/1487)  53.1 W47-W58 N12-N20
      5 1074 1431 (N23W07 L= 28)   256    109        0 (1025/1407-1123/1451) -31.0 W04-W10 N21-N24
      6  711 1376 (N23E17 L=  4)  1032    342        0 ( 584/1300- 868/1431) -30.1 E26-E07 N16-N27
      7  356 1411 (N29E46 L=335)   368     77        0 ( 324/1392- 400/1425)   8.1 E49-E41 N28-N30
      8 1543 1390 (N16W38 L= 59)   328     81        0 (1508/1371-1581/1408) -21.8 W35-W40 N14-N17
      9 1331 1273 (N11W21 L= 42)   200     57        0 (1306/1238-1356/1299) -56.3 W19-W22 N08-N12
     10 1272 1178 (N06W16 L= 37)   505    202        1 (1207/1135-1371/1234) -26.6 W12-W22 N02-N09
     11 1014 1122 (N05W00 L= 21)   264     92        0 (1001/1080-1029/1179) -90.0 W00-W01 N02-N08
     12 1841  940 (S13W58 L= 79)   512     64        0 (1804/ 923-1868/ 965)  15.9 W53-W61 S14-S11
     13 1736  828 (S19W48 L= 69)   616     75        2 (1706/ 797-1771/ 867) -38.7 W44-W51 S21-S16
     14 1662  733 (S24W42 L= 63)   843    230        3 (1640/ 630-1687/ 839) -85.9 W39-W45 S31-S17
     15  188  766 (S08E65 L=315)   596    111        1 ( 170/ 722- 210/ 815)  39.8 E67-E60 S11-S05
     16  866  641 (S23E14 L=  7)   870    188        1 ( 827/ 553- 896/ 725)  61.4 E18-E11 S28-S17
     17  761  512 (S30E24 L=357)   304     61        0 ( 735/ 492- 792/ 534)  14.0 E25-E21 S31-S28
  

    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