Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  759 1649 (N52E10 L=249)   401     72        1 ( 731/1612- 793/1669)  31.0 E12-E06 N49-N54
      2 1599 1545 (N26W55 L=314)   359     72        0 (1572/1517-1628/1569) -63.4 W53-W56 N24-N28
      3 1745 1328 (N10W57 L=316)  1138    194        2 (1687/1270-1809/1408) -65.6 W54-W63 N05-N16
      4  206 1087 (N20E62 L=197)  1264    208        2 ( 142/1050- 292/1153)   2.7 E72-E51 N17-N23
      5 1162  892 (S04W06 L=265)   349    101        2 (1109/ 872-1211/ 905)  -9.5 W03-W09 S05-S02
      6 1294  505 (S31W09 L=268)  3346    652        8 (1019/ 364-1485/ 700)  21.3 E12-W24 S37-S21
      7  517  541 (S17E43 L=216)   297     65        1 ( 499/ 504- 539/ 565)  36.9 E45-E40 S18-S15
      8  740  397 (S29E32 L=227)   505    131        0 ( 723/ 345- 767/ 464)  71.6 E34-E29 S33-S24
      9  532  357 (S30E53 L=206)   789    173        1 ( 507/ 293- 574/ 440)  30.1 E63-E43 S35-S24
     10  809  267 (S41E35 L=224)   423     76        1 ( 776/ 245- 841/ 297)   0.0 E37-E30 S43-S38
  

    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