Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  761 1724 (N33E39 L= 93)   462     87        0 ( 731/1693- 782/1766) -71.6 E41-E37 N30-N36
      2 1212 1673 (N42E06 L=127)  1121    203        1 (1166/1593-1266/1745) -36.9 E12-W00 N37-N47
      3  508 1491 (N12E45 L= 88)  1092    338        4 ( 475/1376- 559/1600) -57.3 E52-E42 N04-N18
      4  850 1545 (N23E25 L=108)   269     61        0 ( 831/1523- 865/1579) -63.4 E25-E22 N21-N25
      5  587 1496 (N14E40 L= 93)   252     60        1 ( 572/1466- 614/1517)  71.6 E41-E38 N12-N15
      6  774 1436 (N15E25 L=107)   317     67        1 ( 759/1396- 792/1472) -59.0 E26-E23 N12-N17
      7 1097 1301 (N14E02 L=130)   221     51        2 (1068/1284-1123/1319)  33.7 E03-W00 N12-N15
      8  186 1247 (S08E67 L= 66)  1561    290        7 ( 135/1197- 223/1296) -56.3 E76-E62 S10-S04
      9  979  843 (S14W02 L=134)   217     43        1 ( 950/ 821- 998/ 861)   0.0 W00-W03 S15-S13
     10  643  640 (S35E15 L=117)   276     40        0 ( 623/ 617- 670/ 668) -45.0 E17-E13 S36-S33
     11 1258  313 (S40W41 L=173)  1278    286        2 (1153/ 285-1398/ 342)  14.7 W31-W51 S42-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