Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  661 1627 (N25E41 L=328)   867    119        0 ( 607/1581- 714/1653) -10.3 E46-E35 N23-N26
      2  511 1496 (N14E46 L=324)   322     30        0 ( 489/1474- 538/1514)  45.0 E47-E43 N12-N15
      3 1467 1313 (N25W21 L= 30)   232     64        0 (1437/1286-1490/1330) -14.0 W18-W23 N23-N25
      4 1578 1030 (N11W34 L= 43)   250     78        0 (1568/ 989-1588/1073) -63.4 W32-W35 N08-N12
      5 1012  761 (S18W05 L= 15)   268     76        0 ( 974/ 739-1054/ 780)  45.0 W03-W07 S20-S16
      6 1256  535 (S26W28 L= 37)   323     91        0 (1211/ 521-1306/ 548)  15.9 W24-W31 S26-S24
      7 1258  395 (S34W35 L= 45)   965    173        2 (1183/ 356-1328/ 444)  -3.8 W27-W42 S35-S32
      8  918  186 (S63W30 L= 39)   488     49        0 ( 907/ 155- 929/ 217) -33.7 W24-W37 S67-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