Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  471 1654 (N19E57 L=169)   407     57        2 ( 445/1635- 500/1677)  26.6 E58-E53 N17-N20
      2  607 1611 (N18E42 L=184)   478    119        1 ( 582/1568- 638/1676) -35.0 E48-E38 N15-N22
      3  315 1519 (N08E63 L=162)  1628    159        5 ( 260/1454- 351/1593) -48.0 E70-E56 N03-N13
      4  503 1489 (N09E44 L=182)   591    143        2 ( 476/1426- 544/1556) -37.9 E48-E39 N05-N12
      5  386 1263 (S07E46 L=180)   356     76        1 ( 344/1242- 422/1284)   0.0 E49-E42 S07-S05
      6  329 1208 (S11E49 L=176)   345     86        0 ( 297/1173- 366/1230) -15.9 E52-E45 S12-S09
      7  905 1159 (S03E10 L=216)   211     38        0 ( 877/1137- 928/1177)   0.0 E11-E08 S03-S01
      8  632 1072 (S13E24 L=202)   206     61        1 ( 620/1035- 643/1106) -76.0 E25-E23 S15-S11
      9  374 1048 (S19E42 L=183)   244     47        1 ( 356/1027- 396/1071)  71.6 E44-E40 S20-S17
     10  477  949 (S23E33 L=193)   395    113        0 ( 443/ 913- 525/ 999) -29.7 E36-E28 S24-S20
     11  225  889 (S31E57 L=169)   475     63        1 ( 203/ 857- 256/ 916) -63.4 E59-E53 S33-S29
     12  206  787 (S38E62 L=164)   499     63        0 ( 173/ 773- 238/ 804)   5.2 E67-E56 S38-S36
  

    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