Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1042 1613 (N39W21 L=161)   909    185        1 ( 997/1529-1093/1685)  33.7 W16-W28 N33-N41
      2 1463 1581 (N23W46 L=185)  1272    179        2 (1388/1530-1537/1614) -41.2 W40-W49 N18-N26
      3  498 1237 (N31E28 L=112)   970    359        0 ( 364/1205- 659/1274) -24.4 E38-E16 N26-N36
      4 1506 1144 (S01W31 L=171)  1066    223        3 (1407/1109-1597/1179) -21.0 W23-W36 S05-N01
      5 1443 1067 (S04W25 L=164)   202     65        0 (1417/1034-1473/1088)  14.0 W22-W26 S04-S02
      6 1660  968 (S16W37 L=176)   657    150        0 (1591/ 937-1706/1000)   0.0 W31-W41 S17-S14
      7 1391  828 (S15W16 L=155)   587    171        0 (1372/ 754-1415/ 927)  60.9 W13-W19 S20-S10
      8 1063  678 (S15E07 L=133)   699    175        0 ( 968/ 640-1145/ 703) -45.0 E12-E03 S19-S11
      9 1165  568 (S24E05 L=135)   366     74        0 (1150/ 530-1180/ 612)  71.6 E06-E03 S27-S21
     10  733  493 (S17E33 L=107)   396    111        0 ( 685/ 471- 795/ 514) -15.9 E36-E28 S18-S15
     11  642  429 (S19E42 L= 98)   419    135        1 ( 611/ 404- 670/ 464) -76.0 E44-E39 S20-S16
     12  770  396 (S25E35 L=105)   489    110        2 ( 737/ 349- 805/ 436)  81.9 E36-E32 S28-S21
     13 1019  212 (S46E32 L=108)  1198    192        3 ( 982/ 157-1053/ 261)  31.6 E38-E25 S50-S41
  

    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