Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  895 1670 (N48E11 L=336)   261     28        0 ( 869/1656- 918/1688)  14.0 E13-E08 N46-N49
      2 1044 1613 (N43W02 L=349)   702    142        1 ( 969/1571-1106/1651) -20.0 E04-W07 N39-N46
      3 1643 1416 (N27W49 L= 35)   811     86        1 (1613/1363-1667/1462)  80.5 W45-W51 N23-N29
      4  783 1273 (N19E16 L=331)   932    283        1 ( 712/1166- 879/1365) -47.3 E21-E09 N12-N25
      5 1426 1168 (N12W26 L= 13)  1022    243        0 (1350/1096-1516/1251) -42.3 W21-W32 N07-N17
      6  488 1160 (N12E36 L=311)   567    101        2 ( 462/1107- 516/1196) -56.3 E38-E33 N08-N14
      7 1210 1114 (N09W12 L=358)   321    128        1 (1185/1052-1230/1167) -74.1 W10-W12 N05-N12
      8 1382  856 (S07W22 L=  9)   242     90        0 (1351/ 833-1419/ 877)  11.3 W20-W25 S08-S05
      9 1285  784 (S12W16 L=  3)   260     70        2 (1255/ 743-1306/ 805) -53.1 W14-W17 S14-S10
  

    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