Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1725 1439 (N21W61 L=116)   470     94        0 (1680/1427-1759/1463)   0.0 W55-W67 N19-N21
      2  132 1020 (N11E74 L=341)  1344    105        3 ( 105/ 978- 155/1060)  39.8 E81-E69 N08-N13
      3 1753  855 (S15W50 L=104)   915    106        2 (1721/ 800-1781/ 904)  45.0 W46-W52 S18-S11
      4 1010  733 (S12E04 L= 51)   223     64        0 ( 987/ 700-1032/ 767)  63.4 E05-E02 S13-S09
      5 1507  688 (S21W28 L= 83)   261     42        0 (1491/ 660-1523/ 715)  63.4 W27-W29 S23-S19
      6 1400  618 (S25W20 L= 75)   295     45        0 (1378/ 591-1423/ 639)  45.0 W18-W22 S26-S23
      7  854  530 (S24E17 L= 38)   293     96        0 ( 803/ 521- 907/ 541)  -8.1 E20-E13 S24-S22
      8  766  464 (S28E25 L= 30)   283    130        1 ( 736/ 441- 811/ 510) -68.2 E27-E22 S29-S24
  

    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