Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  762 1704 (N33E31 L=254)   897    202        0 ( 715/1608- 789/1780) -42.9 E40-E26 N26-N39
      2  468 1633 (N25E52 L=234)   456     65        3 ( 435/1618- 499/1651)   0.0 E56-E48 N23-N26
      3 1219 1612 (N34W04 L=289)   255     85        0 (1194/1573-1242/1646) -51.3 W01-W05 N31-N36
      4  359 1529 (N17E57 L=229)   946    136        2 ( 326/1470- 379/1576) -76.0 E62-E53 N12-N20
      5 1920 1101 (N15W67 L=353)   666     65        0 (1895/1064-1937/1130) -31.0 W63-W70 N13-N16
      6  152  961 (S19E67 L=219)  2870    270        4 ( 111/ 847- 193/1076)  66.8 E75-E61 S25-S11
      7  849  839 (S20E08 L=278)   657    304        2 ( 774/ 772- 949/ 881)  16.7 E12-E01 S22-S17
      8  503  810 (S26E32 L=254)  1031    244        1 ( 407/ 769- 623/ 859)  20.6 E39-E22 S28-S21
      9 1529  602 (S22W42 L=328)  1496    428        1 (1372/ 491-1672/ 742)  58.2 W34-W48 S32-S11
     10  332  708 (S34E49 L=237)   403     37        0 ( 312/ 683- 350/ 737) -90.0 E50-E46 S35-S31
     11 1535  396 (S34W55 L=341)   330     52        0 (1511/ 371-1563/ 420)  76.0 W53-W56 S36-S32
     12  728  288 (S61E14 L=272)   884    134        3 ( 656/ 252- 786/ 316)  -3.2 E24-E06 S64-S58
  

    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