Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  647 1606 (N51E10 L=264)   428     81        1 ( 627/1565- 662/1650)  81.9 E11-E07 N47-N54
      2 1313 1590 (N29W36 L=311)   366     56        1 (1280/1570-1348/1607) -21.8 W34-W39 N26-N29
      3  401 1164 (N28E37 L=237)   447     79        0 ( 359/1127- 432/1202) -50.2 E40-E34 N24-N30
      4  483 1075 (N20E32 L=243)   384     87        0 ( 438/1046- 532/1095)   0.0 E35-E28 N19-N21
      5 1856  789 (S33W56 L=331)  1117    111        0 (1825/ 736-1877/ 852)  66.8 W52-W59 S36-S29
      6 1543  642 (S32W23 L=298)   421    121        0 (1496/ 605-1583/ 692)  20.6 W18-W26 S33-S29
      7 1700  641 (S38W37 L=312)   348     57        0 (1679/ 622-1727/ 669)  71.6 W35-W40 S39-S35
      8 1191  546 (S27E03 L=271)   335    138        0 (1138/ 497-1242/ 585)  14.0 E07-W01 S28-S25
      9 1333  221 (S60E06 L=269)   623     63        2 (1299/ 201-1372/ 244) -35.5 E10-E01 S62-S57
  

    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