Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  953 1833 (N43E38 L=184)   505     35        1 ( 932/1804- 974/1860) -39.8 E41-E35 N40-N45
      2 1329 1693 (N45E03 L=219)  1322    174        2 (1282/1629-1381/1743)  81.3 E04-W01 N38-N52
      3  853 1583 (N22E27 L=194)   252     63        0 ( 815/1572- 890/1597)  21.8 E29-E24 N20-N22
      4 1788 1378 (N41W47 L=269)   899    113        2 (1755/1327-1827/1425) -32.0 W42-W52 N38-N43
      5 1674 1344 (N33W33 L=255)   466     67        0 (1651/1301-1693/1376) -53.1 W31-W35 N30-N35
      6 1471 1342 (N25W17 L=239)   305     94        1 (1425/1314-1514/1364)   9.5 W14-W20 N24-N26
      7  910  701 (S27W03 L=225)   381    142        1 ( 856/ 658- 977/ 735)  54.5 W00-W06 S30-S23
      8 1409  615 (S16W36 L=258)   422    147        0 (1353/ 583-1443/ 677)  74.1 W33-W37 S19-S12
      9 1087  407 (S39W27 L=249)   696    119        1 (1041/ 370-1131/ 454)  90.0 W24-W28 S42-S34
     10 1192  284 (S43W44 L=266)   566    101        1 (1140/ 264-1238/ 308)  71.6 W41-W46 S45-S39
  

    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