Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1458 1576 (N37W39 L= 30)   302     60        1 (1440/1542-1471/1606)  53.1 W37-W41 N34-N38
      2 1396 1477 (N31W30 L= 21)   375     74        2 (1348/1455-1428/1502)   0.0 W26-W32 N28-N32
      3 1550 1358 (N22W39 L= 30)   277     66        0 (1539/1322-1565/1397)  78.7 W37-W40 N19-N24
      4  402 1004 (N05E42 L=310)   940    208        2 ( 371/ 919- 430/1108) -84.8 E44-E39 N00-N11
      5  125 1003 (N04E73 L=278)  3126    186        3 (  97/ 919- 177/1080) -48.8 E81-E64 N00-N08
      6  884  949 (N01E09 L=342)   419    168        0 ( 852/ 883- 900/1033)  69.4 E11-E07 S03-N05
      7  219  789 (S08E62 L=289)  1304    175        2 ( 140/ 736- 284/ 815) -12.5 E74-E54 S11-S06
      8 1047  717 (S15E00 L=351)   346    101        1 ( 992/ 682-1086/ 747) -31.0 E03-W02 S16-S12
      9  333  690 (S14E52 L=300)   317     30        2 ( 314/ 668- 352/ 717) -90.0 E53-E49 S15-S12
     10  336  597 (S21E55 L=296)  1087    178        1 ( 264/ 556- 416/ 647) -19.7 E62-E48 S23-S17
  

    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