Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  794 1727 (N59E03 L=  2)  2095    279        6 ( 696/1651- 911/1789)  10.9 E15-W13 N53-N62
      2 1240 1719 (N46W38 L= 43)  1166    171        3 (1150/1688-1326/1748)  -9.0 W28-W47 N43-N48
      3 1466 1690 (N38W55 L= 61)   802    107        1 (1425/1664-1522/1706) -26.6 W50-W61 N35-N39
      4 1619 1441 (N19W51 L= 56)   323    122        0 (1594/1383-1635/1483)  50.2 W47-W52 N15-N21
      5  645 1341 (N32E20 L=345)   266     82        0 ( 614/1311- 685/1363)  21.8 E23-E17 N30-N33
      6 1374 1262 (N14W26 L= 31)  1213    242        4 (1321/1155-1410/1359)  71.6 W22-W28 N08-N20
      7  630 1090 (N17E24 L=342)   426    135        0 ( 569/1062- 694/1119)   6.3 E28-E19 N15-N17
      8  508  647 (S08E41 L=325)   290     73        1 ( 470/ 627- 546/ 664) -18.4 E43-E37 S09-S07
      9  969  226 (S48E26 L=339)   671     91        1 ( 929/ 204-1016/ 249)  -6.3 E29-E20 S50-S44
  

    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