Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1413 1708 (N33W51 L= 20)   389     60        1 (1385/1693-1446/1721) -11.3 W48-W53 N31-N34
      2 1551 1573 (N21W52 L= 21)   413    107        1 (1506/1555-1584/1591) -15.9 W47-W54 N19-N22
      3 1424 1448 (N19W37 L=  5)   508    133        3 (1380/1404-1466/1505)  18.4 W32-W41 N17-N20
      4 1757 1412 (N06W61 L= 30)   592    140        0 (1718/1398-1804/1427) -16.7 W57-W67 N04-N07
      5 1150 1299 (N19W15 L=343)   259     74        1 (1122/1260-1184/1337) -80.5 W13-W15 N15-N21
      6  541 1006 (N17E30 L=299)   297     72        0 ( 505/ 977- 574/1030) -53.1 E31-E27 N15-N19
      7  997  819 (S04E07 L=322)   190     47        1 ( 977/ 789-1014/ 848) -90.0 E07-E05 S05-S01
      8 1370  622 (S26W10 L=339)   395     82        0 (1329/ 594-1404/ 660)  21.8 W07-W12 S27-S24
      9 1146  175 (S54E26 L=303)  3006    407        6 ( 996/ 105-1324/ 236)  -1.5 E46-E06 S57-S50
  

    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