Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  839 1618 (N49W02 L= 98)   710     94        4 ( 802/1566- 882/1653)  81.9 E01-W05 N44-N51
      2  916 1611 (N46W08 L=105)   255     42        0 ( 896/1584- 933/1631) -90.0 W06-W09 N43-N47
      3  871 1344 (N29E03 L= 94)   526    233        2 ( 809/1273- 962/1398)  20.0 E08-W03 N25-N31
      4  464 1152 (N25E35 L= 62)   787    212        0 ( 376/1087- 519/1210) -39.3 E42-E31 N20-N29
      5  263 1104 (N24E54 L= 42)  1283    284        2 ( 196/1002- 337/1206)  34.5 E63-E46 N18-N30
      6  614 1082 (N18E24 L= 72)   206     96        0 ( 569/1065- 656/1095) -18.4 E27-E21 N16-N19
      7 1109  986 (N03W04 L=101)   186     89        0 (1069/ 962-1140/1021)  21.8 W01-W06 N02-N04
      8  713  661 (S09E26 L= 70)   274     56        1 ( 689/ 640- 740/ 685) -45.0 E27-E24 S10-S06
      9  652  366 (S27E43 L= 54)  1317    226        4 ( 596/ 272- 727/ 458)  71.6 E45-E39 S35-S20
     10  521  429 (S21E50 L= 47)   571    109        0 ( 492/ 397- 567/ 455) -68.2 E53-E47 S23-S18
     11  745  250 (S38E45 L= 52)   399     55        1 ( 721/ 224- 762/ 275)  26.6 E49-E41 S40-S36
  

    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