Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  615 1744 (N30E52 L=355)  1556    333        4 ( 570/1670- 658/1810) -25.2 E60-E43 N25-N34
      2  530 1539 (N15E46 L=  1)  1996    531        6 ( 438/1336- 583/1662) -49.9 E59-E38 N02-N21
      3  378 1539 (N12E59 L=347)  1288    251        4 ( 345/1467- 429/1615) -78.7 E62-E54 N06-N16
      4 1097 1506 (N26E08 L= 38)   946    249        1 (1001/1417-1212/1596) -16.4 E17-W00 N23-N30
      5 1491 1490 (N37W19 L= 65)   468     87        1 (1456/1454-1533/1527) -15.9 W15-W22 N35-N39
      6  674 1434 (N12E32 L= 15)   321    107        0 ( 653/1384- 693/1478) -33.7 E34-E28 N09-N13
      7 1647 1345 (N32W33 L= 80)   299     74        0 (1610/1322-1677/1362)   0.0 W30-W36 N31-N33
      8  577 1235 (S02E32 L= 15)   366     83        2 ( 543/1211- 625/1259)   0.0 E34-E27 S03-N00
      9  725  562 (S39E07 L= 39)   436     85        0 ( 680/ 535- 778/ 587)   7.1 E10-E02 S40-S37
  

    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