Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1124 1877 (N60E25 L= 38)  1691    155        6 (1077/1821-1175/1924) -17.5 E36-E15 N54-N66
      2 1179 1734 (N47E09 L= 55)  2465    297        9 (1088/1608-1258/1815) -33.0 E20-W00 N37-N55
      3  393 1536 (N15E57 L=  6)  2155    269        4 ( 342/1432- 446/1651) -51.3 E68-E53 N07-N22
      4  620 1208 (N00E28 L= 35)  1312    281        2 ( 545/1102- 703/1316)  70.3 E30-E24 S07-N07
      5 1380 1263 (N19W16 L= 80)   226     91        0 (1335/1241-1417/1294)  -8.1 W12-W19 N18-N20
      6  944 1095 (N00E06 L= 57)   236     99        0 ( 928/1045- 954/1145) -63.4 E07-E05 S03-N02
      7 1196 1008 (N00W10 L= 74)   277    148        0 (1182/ 935-1220/1069) -74.1 W09-W12 S03-N04
      8  849 1002 (S07E10 L= 54)   261    102        0 ( 807/ 968- 889/1029)  51.3 E11-E07 S09-S04
      9 1792  770 (N00W59 L=122)  3619    305        5 (1727/ 635-1849/ 918) -31.8 W49-W70 S06-N07
     10 1452  855 (S03W29 L= 92)   543    103        1 (1413/ 822-1495/ 893)  59.0 W26-W31 S05-N00
     11 1182  171 (S54W52 L=115)   507     67        0 (1157/ 143-1204/ 201) -31.0 W47-W58 S56-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