Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  832 1616 (N49W01 L= 98)   741    103        2 ( 792/1562- 867/1650) -90.0 E02-W04 N44-N52
      2  841 1326 (N29E05 L= 92)   407    148        2 ( 807/1272- 896/1383)  32.0 E08-W00 N25-N31
      3  439 1149 (N25E37 L= 60)  1583    524        2 ( 294/1065- 643/1226) -23.2 E51-E22 N16-N31
      4  228 1067 (N22E59 L= 38)   842    145        1 ( 194/1002- 263/1126)  37.9 E63-E54 N18-N25
      5  706  660 (S09E27 L= 70)   233     55        1 ( 684/ 638- 736/ 682) -56.3 E28-E25 S10-S06
      6 1056  477 (S27E10 L= 87)   271     54        1 (1032/ 449-1080/ 506)  36.9 E12-E08 S29-S25
      7 1283  436 (S36W04 L=101)   276     54        2 (1257/ 407-1308/ 469) -59.0 W02-W05 S38-S33
      8  511  431 (S21E51 L= 46)   517     63        0 ( 489/ 401- 542/ 455)  63.4 E53-E47 S22-S18
      9  646  366 (S27E43 L= 54)  1327    220        3 ( 595/ 276- 722/ 454)  75.1 E45-E40 S35-S20
     10  741  245 (S39E46 L= 51)   402     55        1 ( 722/ 218- 756/ 274)  32.0 E50-E42 S41-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