Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1379 1799 (N39W56 L=151)  1340    176        2 (1312/1758-1422/1829) -18.4 W50-W61 N36-N41
      2 1339 1492 (N24W34 L=129)   223     50        1 (1307/1478-1366/1506) -33.7 W31-W35 N22-N24
      3  765 1304 (N29E08 L= 87)   531    142        1 ( 715/1262- 810/1362) -69.4 E10-E04 N26-N34
      4 1762 1086 (S11W48 L=143)  3730    569        6 (1714/ 855-1808/1328)  58.0 W40-W55 S23-N02
      5 1126 1131 (N10W08 L=104)   345    142        0 (1049/1090-1191/1159) -40.6 W04-W11 N06-N12
      6  704 1051 (N16E18 L= 77)   273     90        0 ( 667/1032- 753/1077) -38.7 E19-E14 N14-N18
      7  235 1002 (N23E54 L= 41)  1688    348        3 ( 143/ 954- 328/1035) -15.5 E64-E46 N19-N26
      8  762  513 (S16E30 L= 65)  1166    310        2 ( 726/ 387- 807/ 634)  63.4 E35-E23 S24-S08
      9  540  393 (S18E51 L= 44)   427     97        1 ( 500/ 363- 581/ 433) -71.6 E53-E48 S21-S15
     10  913  374 (S29E27 L= 68)   492    114        2 ( 852/ 361- 975/ 394) -23.2 E30-E23 S30-S26
  

    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