Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  969 1800 (N40E32 L=258)   316     39        0 ( 949/1780- 988/1823) -63.4 E33-E30 N37-N41
      2  736 1733 (N28E42 L=247)   563    172        0 ( 707/1680- 761/1784) -31.0 E47-E37 N24-N30
      3  591 1473 (N08E39 L=251)   265     69        1 ( 569/1442- 610/1505) -21.8 E41-E36 N05-N09
      4 1064 1326 (N11E06 L=283)  1030    322        0 (1035/1159-1087/1466) -59.5 E11-E01 N01-N19
      5 1282 1070 (N03W13 L=302)   363    128        0 (1244/1019-1306/1129) -71.6 W10-W14 N00-N06
      6  228  991 (S28E54 L=236)  1881    398        0 ( 129/ 936- 331/1074)  22.8 E64-E45 S32-S21
      7  753 1006 (S15E15 L=275)  1310    348        1 ( 572/ 957- 906/1039)  15.9 E26-E05 S17-S11
      8 1423  342 (S31W52 L=342)   414     35        1 (1401/ 321-1444/ 364) -71.6 W49-W53 S32-S29
      9 1277  313 (S38W45 L=334)   993    165        2 (1200/ 285-1344/ 337)  15.3 W38-W49 S40-S35
  

    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