Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  819 1863 (N40E52 L=348)  1156    144        0 ( 780/1794- 843/1908) -23.2 E64-E43 N34-N44
      2  401 1628 (N13E65 L=334)  1982    262        2 ( 352/1532- 462/1751) -54.8 E74-E59 N06-N23
      3 1169 1611 (N32E10 L= 29)   446    103        2 (1147/1568-1191/1665) -49.4 E13-E07 N29-N36
      4  428 1536 (N08E55 L=345)  2410    293        3 ( 315/1454- 545/1631)  49.4 E63-E50 N01-N15
      5 1113 1296 (N11E03 L= 37)   242    158        1 (1072/1240-1160/1354)  77.5 E04-E01 N06-N15
      6  218 1167 (S18E61 L=339)   955     89        2 ( 175/1130- 256/1205)  24.0 E66-E56 S20-S15
      7  532  877 (S28E27 L= 13)   691    195        1 ( 454/ 827- 616/ 943)  47.7 E32-E21 S32-S21
      8 1325  621 (S19W30 L= 70)  1474    483        5 (1145/ 469-1460/ 740)  76.5 W26-W34 S33-S08
      9 1056  689 (S24W12 L= 51)   207    119        0 ( 996/ 682-1115/ 700)  29.7 W08-W15 S25-S21
     10  834  654 (S32E00 L= 39)   452    160        1 ( 760/ 633- 915/ 675)  36.9 E04-W04 S35-S28
     11  907  505 (S39W10 L= 50)  1227    349        1 ( 788/ 426-1056/ 584)  59.9 W04-W16 S48-S29
  

    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