Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  651 1608 (N51E10 L=264)   385     67        1 ( 628/1569- 667/1640)  59.0 E11-E07 N47-N53
      2 1323 1593 (N28W37 L=311)   351     82        0 (1285/1575-1369/1606) -18.4 W34-W40 N27-N29
      3  400 1172 (N28E37 L=237)   462    110        1 ( 346/1133- 437/1206) -35.5 E41-E34 N25-N31
      4  491 1078 (N20E31 L=243)   387     93        0 ( 445/1046- 541/1098)   7.1 E35-E27 N19-N21
      5  202  928 (N18E61 L=213)   406     68        2 ( 183/ 895- 227/ 953)  36.9 E62-E56 N15-N19
      6 1860  787 (S34W57 L=331)  1241    117        0 (1834/ 735-1882/ 852)  66.8 W53-W60 S37-S29
      7 1702  642 (S37W38 L=312)   408     50        1 (1681/ 621-1731/ 673)  71.6 W35-W41 S39-S35
      8 1524  619 (S33W21 L=295)   242     72        0 (1485/ 604-1555/ 639)   0.0 W17-W23 S33-S31
      9 1192  544 (S27E03 L=271)   309    127        0 (1142/ 497-1238/ 583)  23.2 E07-W00 S29-S25
     10 1338  223 (S60E05 L=269)   596     62        1 (1304/ 206-1374/ 246) -39.8 E08-E01 S63-S56
  

    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