Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1538 1684 (N34W57 L= 55)   410     85        2 (1519/1653-1554/1711)  53.1 W54-W59 N31-N35
      2 1722 1365 (N11W56 L= 54)   997    161        2 (1696/1293-1755/1443)  78.7 W53-W58 N05-N15
      3 1560 1311 (N12W40 L= 39)   321     46        1 (1533/1286-1587/1334) -56.3 W38-W42 N10-N13
      4 1366 1239 (N12W25 L= 23)   219     86        0 (1321/1228-1412/1258) -38.7 W22-W28 N10-N14
      5  199 1121 (N25E59 L=299)  1412    234        3 ( 169/1033- 230/1202) -68.2 E62-E56 N19-N29
      6  782 1102 (N16E13 L=345)   549    183        1 ( 727/1018- 849/1152)  29.1 E17-E08 N11-N18
      7  575 1101 (N19E27 L=331)   683    179        2 ( 489/1056- 643/1141) -31.0 E32-E22 N15-N22
      8  916 1085 (N12E05 L=353)   234    108        0 ( 858/1050- 963/1131) -60.3 E08-E03 N09-N16
      9  719  862 (N03E21 L=337)   337     96        1 ( 678/ 832- 763/ 885)   9.5 E24-E18 N01-N03
     10  650  650 (S09E30 L=328)   765    237        0 ( 577/ 603- 764/ 722) -35.0 E33-E22 S12-S03
     11  341  487 (S16E63 L=295)  1583    293        2 ( 314/ 386- 367/ 577)  65.0 E67-E58 S24-S09
  

    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