Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1779 1450 (N08W64 L=310)  2171    348        7 (1708/1359-1849/1535) -57.5 W57-W71 N01-N12
      2 1773 1377 (N04W60 L=306)   383     79        2 (1749/1350-1801/1413) -76.0 W58-W62 N02-N06
      3  133 1174 (N31E62 L=183)  1118    212        2 ( 100/1124- 162/1213) -50.2 E66-E58 N28-N34
      4  731 1023 (N10E17 L=229)   289     79        1 ( 684/ 992- 758/1047)   0.0 E19-E14 N08-N11
      5  857  908 (N01E12 L=234)   180     74        0 ( 822/ 882- 899/ 938) -53.1 E13-E10 S01-N03
      6  942  910 (S01E07 L=239)   334     74        0 ( 899/ 889- 979/ 930)   0.0 E09-E04 S02-N00
      7  221  627 (S02E69 L=176)  2549    355        2 ( 156/ 541- 307/ 700) -42.0 E81-E62 S08-N02
      8  330  448 (S15E67 L=179)  2009    299        3 ( 299/ 381- 362/ 498)  48.8 E73-E61 S19-S11
      9  562  442 (S20E45 L=201)   332     54        0 ( 534/ 422- 587/ 457) -33.7 E46-E42 S21-S18
  

    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