Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  576 1666 (N19E49 L=225)   502     82        1 ( 542/1638- 616/1698)  80.5 E51-E47 N16-N22
      2  909 1665 (N28E27 L=247)   290     87        0 ( 878/1625- 935/1698) -33.7 E30-E24 N25-N29
      3 1218 1632 (N35E07 L=267)   263     45        1 (1203/1607-1239/1655) -45.0 E09-E05 N33-N36
      4  354 1414 (S01E55 L=219)   467     65        0 ( 338/1384- 368/1448) -56.3 E56-E52 S02-N01
      5  827 1373 (N08E21 L=253)   588    208        0 ( 756/1311- 919/1445) -12.1 E27-E13 N06-N10
      6 1291 1226 (N12W09 L=283)  1170    291        1 (1254/1058-1314/1361) -58.0 W03-W14 N03-N19
      7  201 1196 (S16E62 L=213)   832    136        2 ( 171/1163- 238/1232)  26.6 E65-E57 S18-S13
      8 1498  947 (N03W29 L=303)   389    117        0 (1463/ 896-1518/1010) -80.5 W26-W30 N00-N06
      9  341  924 (S30E42 L=232)  1014    296        0 ( 241/ 891- 438/ 970)  16.4 E51-E34 S32-S24
     10  993  886 (S15W02 L=276)  1110    336        0 ( 817/ 847-1153/ 922)  17.5 E08-W11 S17-S11
     11  208  682 (S48E56 L=218)   553     80        0 ( 187/ 667- 239/ 705)   0.0 E58-E52 S49-S46
     12 1432  349 (S30W52 L=326)   375    101        0 (1376/ 324-1469/ 373)  63.4 W49-W54 S33-S27
  

    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