Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  530 1663 (N63E27 L=276)   364     87        0 ( 487/1655- 574/1675)  -7.6 E34-E19 N61-N65
      2 1529 1546 (N22W50 L=353)   327     50        0 (1502/1525-1559/1566) -63.4 W48-W52 N19-N23
      3  429 1330 (N39E35 L=268)   683    118        0 ( 410/1278- 458/1392) -76.0 E37-E32 N35-N43
      4  738 1278 (N29E11 L=292)   525    144        1 ( 692/1247- 808/1316) -40.6 E14-E06 N25-N32
      5 1758  697 (S34W44 L=347)   947    151        1 (1717/ 644-1806/ 763)  29.1 W39-W49 S36-S31
      6 1380  589 (S28W10 L=313)   285     81        0 (1364/ 545-1397/ 627)  80.5 W09-W11 S31-S25
      7  735  355 (S26E39 L=264)   484     93        1 ( 714/ 314- 761/ 397)  63.4 E41-E36 S29-S23
  

    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