Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  960 1831 (N43E38 L=184)   522     45        0 ( 938/1802- 983/1856) -33.7 E40-E34 N40-N45
      2 1336 1693 (N46E02 L=219)  1335    173        3 (1293/1629-1393/1744)  81.9 E04-W01 N38-N52
      3 1794 1374 (N41W48 L=269)   809     96        3 (1764/1331-1830/1410) -23.2 W43-W52 N38-N43
      4 1680 1341 (N33W34 L=256)   497     63        0 (1660/1294-1698/1377) -59.0 W32-W36 N30-N35
      5 1472 1342 (N25W17 L=239)   265     86        0 (1432/1327-1519/1363)   8.1 W14-W21 N24-N26
      6  909  693 (S27W03 L=225)   390    144        0 ( 857/ 648- 974/ 728)  58.0 W01-W06 S31-S23
      7 1413  602 (S17W37 L=258)   353    114        0 (1360/ 578-1439/ 644)  59.0 W33-W38 S19-S14
      8 1092  403 (S39W28 L=249)   737    118        2 (1042/ 365-1141/ 453)  82.9 W25-W29 S42-S34
      9 1233  313 (S39W44 L=265)   919    231        0 (1152/ 257-1321/ 392)  90.0 W41-W47 S45-S31
     10  833  297 (S57W18 L=239)   305     57        0 ( 821/ 266- 850/ 322) -50.2 W15-W21 S59-S53
  

    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