Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  743 1608 (N20E33 L=324)   283     73        2 ( 713/1570- 762/1631) -45.0 E35-E31 N17-N21
      2  674 1564 (N16E36 L=322)   636     90        2 ( 631/1537- 730/1585)  33.7 E39-E32 N14-N18
      3 1412 1370 (N23W14 L= 12)  1066    355        0 (1347/1213-1472/1523) -52.6 W07-W21 N14-N31
      4 1671 1191 (N20W36 L= 34)  1054    295        1 (1594/1069-1741/1285) -29.7 W28-W43 N14-N23
      5  446  895 (S28E35 L=323)   274     46        0 ( 427/ 858- 458/ 922) -56.3 E36-E33 S29-S26
      6 1026  593 (S32W12 L= 10)   249     41        0 ( 997/ 575-1057/ 608)   0.0 W10-W14 S32-S30
      7  896  529 (S39W06 L=  4)   646    182        0 ( 833/ 470- 974/ 560)  56.3 W03-W10 S45-S35
      8  397  446 (S61E47 L=311)   669     62        1 ( 364/ 420- 422/ 480)  16.4 E57-E40 S63-S57
      9  546  385 (S62E23 L=335)   536     61        0 ( 510/ 361- 585/ 412)   0.0 E28-E17 S63-S59
  

    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