Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  892 1777 (N60W11 L=287)   421     58        1 ( 867/1743- 913/1802)  50.2 W07-W14 N56-N62
      2 1439 1599 (N32W45 L=321)  1547    240        3 (1372/1557-1551/1670) -42.0 W38-W53 N27-N37
      3  238 1455 (N45E71 L=205)  1078    105        2 ( 203/1425- 286/1478)   7.1 E83-E61 N42-N46
      4  580 1455 (N41E25 L=251)  1072    270        1 ( 452/1427- 699/1476) -14.6 E38-E15 N37-N43
      5 1724  660 (S31W44 L=321)  2862    408        4 (1601/ 558-1848/ 835)  19.8 W32-W57 S36-S22
      6  208  654 (S05E71 L=205)   880     78        2 ( 184/ 606- 232/ 685)  59.0 E76-E67 S08-S03
      7 1423  633 (S25W18 L=295)   402    106        0 (1361/ 622-1473/ 649) -14.0 W14-W22 S26-S23
      8  621  447 (S23E39 L=237)   463     85        0 ( 585/ 415- 660/ 477) -53.1 E41-E36 S25-S20
      9  998  252 (S45E21 L=255)  3710    498        5 ( 785/ 161-1211/ 346)   8.5 E49-E01 S50-S40
  

    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